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, We use a density functional theoretical approach to calculate the pair distribution function and 

the effective interactions in homogeneous fluids of spinless charged bosons. The scheme involves the 
\ self-consistent solution of a two-particle scattering problem with an effective scattering potential 

. which embodies many-body effects and is adjusted to the compressibility sum rule. Numerical 

! i ' results are presented over an extensive range of density in both three and two dimensions. 

o : 

Exchange and correlations in fluids of charged particles have been a focus of interest in many-body physics for many 
decades. An important manifestation of these effects is the equilibrium distribution g(r) of pairs of particles in the 
ground state of the homogeneous fluid, which describes its state of short-range order. Knowledge of g(r) as a function 
of the fluid density determines its energy and is essential in the construction of energy functionals for applications of 
density functional theory (DFT) beyond the local density approximation 1 . Main attention has been given to electron 
fluids because of their relevance in the physics of metals and semiconductors, but there also is an interest for plasmas 
of charged bosons in quantum statistical mechanics, for instance in regard to condensates of point-like Cooper pairs 
as models for superfluid states 2 or to the equation of state and nuclear reactions in astrophysical plasmas 3 . 

There has recently been a renewed interest in the study of g(r) for electron gas models within a two-body scattering 
approach stemming from work by Overhauser 4 . In brief, g(r) is obtained from the solution of a Schrdinger equation for 
particle-pair wave functions with effective scattering potentials which, starting from a simple electrostatic model 4,5 , 
have been developed into a self-consistent Hartree model 6 and into spin-dependent effective pair interactions . In the 
present work we derive a DFT basis for such an approach and, using an earlier modelling of the effective interactions 
^ . to incorporate the thermodynamic sum rules 8 , we develop a fully self-contained and self-consistent determination of 
g(r) and of the effective scattering potential. Our focus here is on plasmas of spinless Bose particles, leaving aside the 
\q ' further complications that arise in the presence of the spin degree of freedom as already briefly discussed elsewhere 8 ' 9 . 

We consider, therefore, a quantum fluid of point-like bosons having charge e and mass m, which are free to move 
in dimensionality D = 2 or 3 at zero temperature and are neutralized by a uniformly charged background of density 
n. This model will be referred to as a charged Bose fluid (CBF). Our aim is to use DFT for building a self-consistent 
theory that gives the pair distribution function g(r) as output. From the Hohenberg-Kohn theorem 10 the ground 
state energy of the fluid in the presence of an external potential V cx t ( r ) can be written as a functional of the density 
I ■ profile n(r) in the form 
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E gs [n(r)} = T s [n(r)] + E u [n(r)} + / d D r V CKt (r)An(r) + E D [n(r)] W 



where An(r) = n(r) — n and T s is the ideal kinetic energy functional. The Hartree term is given by 



%[»(r)] = - / d D r I d D r'v(\r~r'\) An(r)An(r') (2) 



where v(\r — r'|) = e 2 /|r — r'\. The last term in Eq. (1) is the correlation energy functional, which contains all the 
quantum many-body (QMB) effects. In Eqs. (1) and (2) the presence of a neutralizing background has been taken 
into account. 

The quantity n[g{r) — 1] in the homogeneous fluid can be viewed as the distortion in the density profile that 
instantaneously surrounds a particle of the fluid located at position r = 11 . As usual g(r) is defined by setting 
ng(r)Q,£)r D ^ 1 dr, with = 2tt and = 4ir, equal to the average number of particles inside a shell of radius r and 
thickness dr centered on the particle at the origin. The appropriate ground-state energy functional for the surrounding 
fluid is obtained from Eq. (1) with the formal replacements V^ x t(r) — > v(r) and 
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An(r) = n[g(r) - 1] 



(3) 



Following the treatment given for inhomogeneous fluids in the book of Dreizler and Gross 10 , a formal expression for 
the QMB energy functional can be obtained via an adiabatic connection formula 10 , 



£7 c [n(r)] = ~J dX J d D r J d D r' v(\r - r'|) n(v)n(v') [<?f(r,r') - 1 



(4) 



where the integration over the coupling strength A accounts for the shift in kinetic energy that accompanies the 
switching on of the interactions. In Eq. (4) we have n(r) = ng(r) and gf \r,r') measures the probability of finding 
two particles at r and r' when a third particle is at the origin, the interaction potential being v\(\r — r'|) = Ae 2 /|r — r'|. 
Of course, g x (r, r ) depends functionally on n(r). 

An Euler-Lagrange equation for g(r) in the boson fluid can now be obtained from the variational principle, using 
the ideal kinetic energy functional which for bosons is given by the von Weizsacker expression 12 . This is 
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8m r J ' ' g(r) 

where m r = m/2 is the reduced mass of a particle pair. Taking the zero of energy at the chemical potential, the 
equation for g(r) reads 

^-v 2 + y KS (r)| V^M = o (6) 



where V KS (r) is the Kohn-Sham scattering potential, 

V Jta (r)=v(r) + J d D r'v( 



r'l)An(r') 
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In fact, Eq. (6) can also be obtained from the Kohn-Sham mapping in DFT 10 by building An(r) from Kohn-Sham 
scattering orbitals <I>k(r). Since at zero temperature all bosons in the reference Kohn-Sham ideal gas are in the zero 
momentum state, we have g(r) oc | <&k=o ( r ) 1 2 an d Eq. (6) is just the Kohn-Sham Schrodinger equation for the pair 
wave function $k=o( r ) at zero relative momentum. We conclude, therefore, that for a Bose fluid the scattering-theory 
approach to the pair distribution function admits a rigorous DFT derivation, which yields Eqs. (6) and (7). 

Of course, the functional dependence of the QMB energy on density is not known and at this point we have to 
resort to some approximations. Their goodness can only be gauged a posteriori. Firstly, the three-body correlation 
function gf\r, r') in Eq. (4) would lead us into a hierarchy of higher-order correlation functions. The simplest way 
of truncating this hierarchy is to replace g^\r,r') by g\(\r — r'|), in analogy with what has been done in treating 
the equation of motion for the Wigner distribution function in the electron gas subject to external potentials 13 . This 
approximation can be expected to work well at strong coupling where the probability of simultaneously finding three 
particles inside a radius r s aB is small, such a "strong coupling" situation being reached at lower values of r s in 2D 
than in 3D. As a second approximation we expand the QMB energy in a functional Taylor series in powers of An(r) 
up to second order terms. With the definition 



f(|r 



S 2 E c [n(r 



6n(r)6n(r') 



Ara(r)=0 



we find in Fourier transform 



V KS (q)=v(q)+v(q)[l-G(q)} An(q) 



(8) 



(9) 



Here, G(q) = —f(q)/v(q) is the so-called local field factor, defined in terms of the Fourier transform f(q) of the 
QMB kernel in Eq (8), and v(q) is the Fourier transform of the Coulomb potential (i.e. v(q) = Aire 2 /q 2 in D = 3, 
v(q) = 2ire 2 /q in D = 2). Finally, An(q) in Eq. (9), is given by 
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An(q) = S(q) - 1 



(10) 



where S(q) is the structure factor, related to the pair function by 

S(q) = l + n J d D r [g(r) - 1] exp (-iq • r) . (11) 

The solution of Eqs. (6) and (9) can therefore be carried out self-consistently, given knowledge of the local-field factor 
G(q). The results of such a self-consistent scheme (SCS) for the 3D- CBF, using the data on G(q) from the Quantum 
Monte Carlo (QMC) study of Moroni et al. 14 , will be reported in Figures 1 and 2 below. 

In fact, the present approach can be extended into a fully self-contained theory in which the local field factor is 
self-consistently determined over the relevant g-range during the calculation rather than taken as input from QMC 
data. In such a fully self-consistent scheme (FSCS) we adopt a closure relation between G(q) and S(q), which satisfies 
the compressibility sum rule 8 . That is, we set 

G(q)=D n g(q) (12) 

where the differential operator D n is defined by 
while Q{q) is given by 

m = -^h( dX !w^ ^ + - 1] ' (14) 

with S\{q) being the partial structure factor at coupling constant A. Although these expressions are strictly correct 
only in the long-wavelength limit, they yield a good account of QMC data on local field factors over the relevant range 
of q (see Ref. 8 and Figure 4 below). The explicit proof that Eqs. (12)-(14) satisfy the compressibility sum rule at 
long wavelengths requires two simple steps: (i) the dependence of the integral in Eq. (14) on wave number q can be 
neglected in the limit q — > 0, yielding Q{q) — > —2e c (n)/[nv(q)] in this limit where e c (n) is the correlation energy of 
the CBF; and hence (ii) D n Q(q) — > —l/{n 2 Kv(q)] with k being the compressibility of the CBF. 

We turn at this point to a presentation of our numerical result. We have solved the FSCS based on Eqs. (6), (9)-(ll) 
and (12)-(14) for a 3D-CBF with coupling strength up to r s = 20 (with r s a B = (47m/3) _1/3 ) and for a 2D-CBF with 
coupling strength up to r s = 10 (with r s aB = (7rn) -1 / 2 ). The main results of our work are shown in Figures 1-4. 

In Figure 1 we show that our FSCS results for g(r) in the 3D-CBF at r s = 10 and r s = 20 arc in excellent agreement 
with the QMC data of Moroni et al. 14 . We also compare our results with those that we have obtained in the Hartree 
approximation (HA, setting G(q) — 0) and with those that we have recalculated by the hypernetted-chain Euler- 
Lagrange approach (HNC/EL) of Apaja et al. 15 , the latter being still based on Eq. (6) but with a different choice for 
the dependence of the scattering potential on S(q). The HA is unable to reproduce the emergence of a first-neighbor 
shell with increasing coupling strength, as already noted in its use for the study of pair correlations in the electron 
gas 6 . The two self-consistent approaches to the dependence of the scattering potential on the structure factor, on the 
other hand, are in very good agreement with each other. 

The inset in Figure 1 shows that our SCS results at r s — 20 are in very good agreement with those obtained in the 
FSCS, stressing that our self-consistent determination of G(q) from the compressibility sum rule also accounts with 
sufficient accuracy for this function over the relevant range of wave number. We have made use of this observation 
to test the accuracy of the present theory for g(r) at very strong coupling, as is shown in Figure 2 by reporting SCS 
results for the 3D-CBF at r s = 100 in comparison again with the QMC data of Moroni et al. 14 on g(r) and with the 
results of the HA and of the HNC/EL 15 . 

In Figure 3 we illustrate the quality of our results for the 2D-CBF at r s — 5 and r s = 10. No data are as yet 
available for the local-field factor in this model system, and we show in this Figure our FSCS results in comparison 
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with QMC data on g(r) kindly sent to us by Dr. Moroni prior to publication. Again we obtain excellent agreement 
with the QMC data and with the HNC/EL theory of Apaja et al. 15 . 

Finally, in Figure 4 we show our FSCS results for the local- field factor G(q) in the 3D-CBF at r s = 10, in comparison 
with the QMC data of Moroni et al. 14 , and for the 2D-CBF at r s = 5. Deviations from the QMC data are seen to 
emerge in 3D with increasing wave number, starting at qr s as > 1.5 and becoming very large at qr s aB > 4. However, 
as already stressed in the discussion of the results for g(r) in Figure 1, these discrepancies are scarcely of any relevance 
in the calculation of pair correlations once the compressibility sum rule is embodied into the theory. 

In summary, we have shown that a scattering-theory approach to pair correlations in boson plasmas has a sound 
theoretical justification within the framework of DFT and leads to fully quantitative results when self-consistence 
between pair correlations and effective particle-particle interactions is incorporated into the theory. It may be worth 
exploring in the future alternative approximations to the truncated expansion of the QMB energy functional. 

As a final remark we wish to point out that, having taken the Kohn-Sham viewpoint in our choice of an ideal Bose 
gas as the reference DFT fluid, our approach carries no information on momentum distributions and in particular 
on the depletion of the condensate which occurs with increasing coupling strength. An extended version of the DFT 
scheme, which adopts both the particle density and the order parameter of the condensate as basic variables of the 
inhomogeneous fluid and uses the Bogoliubov - de Gennes equations as reference, has been developed for such wider 
purposes 16 . 
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FIG. 1. The pair distribution function g(r) in a 3D-CBF at r s = 10 and 20, as a function of r in units of r s os. In the main 
body of the Figure the results of the FSCS (full lines), of the HA (dotted lines) and of the HNC/EL theory (dashed lines) are 
compared with QMC data from Ref. 14 (crosses). The curves at r s = 20 have been shifted upwards by 0.4. In the inset the 
SCS results at r s — 20 (dotted line) are compared with the FSCS results (solid line) and with the QMC data (crosses). 
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FIG. 3. The pair distribution function g(r) in a 2D-CBF at r s — 5 and 10, as a function of r in units of r s as. The results 
of the FSCS (full lines), of the HA (dotted lines) and of the HNC/EL theory (dashed lines) are compared with QMC data of 
Dr. S. Moroni (crosses, unpublished). The curves at r s — 10 have been shifted upwards by 0.4. 
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FIG. 4. The local-field factor G(q) of the 3D-CBF at r s = 10 and of the 2D-CBF at r s = 5, as a function of q in units of 
(r s as) _1 . The FSCS results are shown as full lines and are compared for the 3D-CBF with QMC data from Ref. 14 (crosses) 
and with an interpolation formula to the QMC results reported in Ref. 14 (dashed line) . The curve giving the FSCS results for 
the 2D-CBF has been shifted upwards by 0.5. 
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